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ABSTRACT 

We  study  a  natural  clasa  of  appropriate  viscosity  matrices  for 
strictly  hyperbolic  systems  of  conservation  laws  in  one  space  dimension, 
ut  +  f(u)x  »  0,  u  e  iP.  These  matrices  are  admissible  in  the  sense  that 
small  amplitude  shock  wave  solutions  of  the  hyperbolic  system  are  shorn  to  be 
limits  of  smooth  traveling  wave  solutions  of  the  parabolic  system 
u,.  +  f(u)  *■  v(Du  )  as  v  ♦  0  if  D  is  in  this  class.  The  class  is 
determined  by  a  linearised  stability  requirements  The  Cauchy  problem  for  the 
equation  ufc  +  f 1 (uQ )ux  -  vDu^  should  be  well  posed  in  L2  uniformly  in 
v  as  v  >  0.  Previous  examples  of  Inadmissible  viscosity  matrices  are 
accounted  for  through  violation  of  the  stability  criterion. 
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\  SIGNIFICANCE  AND  EXPLANATION 

M 

Many  aquations  of  mathematical  physics  take  the  form  of  nonlinear 
hyperbolic  systems  of  conservation  laws.  With  small  dissipative  effects 
neglected,  typically  smooth  solutions  must  develop  discontinuities  (shocks  in 
finite  time.  Re incorporating  dissipation  helps  select  those  discontinuities 
which  are  physically  meaningful.  For  this  purpose,  many  different  sorts  of 

dissipation  will  do j  in  particular,  the  physical  viscosity  is  typically 

'r<~  /"  J 

degenerate  and  not  convenient.  In  this  paper  we  provide  a  thorough 
understanding  of  what  sorts  of  second  order  viscosity  terms  smooth  the 
physical  discontinuities.  A  natural  class  of  '"admissible'*  viscosity  terms  is 


determined  based  on  a  simple  linearized  stability  condition. 

A 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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|1.  INTRODUCTION 


The  siaplest 

discontinuous  solutions  of 

th« 

n  x  ■  system  of  hyperbolic 

laws 

(1.1) 

“t  *  *<«>x 

-  0, 

u  e  ^ 

are  the  shock  wave 

solutions  defined  by 

X  <  St 

(1.2) 

u(x,t)  -  - 

“r 

X  >  St 

where  and  ara  constant  vectors  which  together  with  the  constant  a  satisfy  the 

RanXine-Hugonlot  conditions 

(1.3)  -stUjj  -  ♦  f<u„)  -  ftu^)  -  0 

and  a  suitable  strict  entropy  condition!  Lax's  shock  inequalities  in  the  genuinely 
nonlinear  case  (see  (1.15)  below),  and  Liu's  strict  condition  (S)  in  the  general  case  (see 
section  3). 

We  assuaa  that  the  systan  (1.1)  is  strictly  hyperbolic.  Thus,  if  A(u)  -  3f/Ju  is 
the  nxn  Jacobian  Matrix,  A(u)  has  «  distinct  real  eigenTalues,  ordered 
1^(0)  <  X2(u)  <  ...  <  X^(u)  with  corresponding  right  and  left  eigenvectors  r^tu)  and 
**(<»)  for  },k  ■  1,...,a,  satisfying 

A<“,rj  *  Vj  »*<“>Sc  *  Sc1)! 

(1.4) 

Vrj  '  6kj 

An  eigenvalue  A^(u)  is  called  genuinely  nonlinear  if  71^  «r^<u)  never  vanishes. 
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In  1959  Gelfand  introduced  the  following  probleai  Show  that  a  discontinuous  solution 
of  the  form  (1.2)  (so  satisfying  (1.3))  is  the  Halt  of  special  saooth  solutions 
uw  _  o( —  j  of  a  reasonable  parabolic  systea 

(1.5)  o'*  ♦  f(uV)  ■  v(o(uV)uU)  as  u  ♦  0 

t  X  XX 

if  (and  only  if)  the  entropy  condition  is  satisfied.  In  particular,  one  should  determine 
the  class  of  viscosity  aatrlces  D(u)  for  which  the  above  is  trues  such  a  matrix  is 
called  acbalaelble.  Use  saooth  traveling  wave  solution  o(x-  ■  )  is  called  a  viscous 

shock  profile. 

The  existence  of  the  traveling  wave  o(— at )  having  the  desired  Halt  (1.2) 

X  "•  ft 

requires  that,  with  5  “  — - - ,  0(C)  should  satisfy  the  a  x  a  systea  of  nonlinear 

ODK's 

(1.6)  D(u)ue  -  -s(0  -  u^)  ♦  f(0)  -  f(uL) 
for  ■•<(<•  together  with  the  boundary  conditions 

(1.7)  lia  0(C)  -  ul  11a  0(C)  -  u„ 

^  R 

That  it,  the  autonomous  systea  of  ODK's  (1.6)  should  admit  a  trajectory  connecting  the 
critical  point  on  the  left  to  the  critical  point  u^  on  the  right. 

Several  authors  have  established  sufficient  conditions  ensuring  the  existence  of  the 
connecting  orbit  for  (1.6),  (1.7)  under  various  asswptions,  typically  including  either 
a  •  2  or  | u^  -  uR|  aaall  (weak  shocks),  and  either  genuine  nonlinearity  or  D(u)  5  I 
(see  (4],  (1],  (2), (3),  [11]).  h  notable  exception  is  Nock's  more  recent  paper  [12],  which 
finds  a  broad  class  of  aatrlces  0(u)  admissible  for  general  a  and  strong  shocks,  asking 
global  assumptions  of  genuine  nonlinearity  and  the  existence  of  a  convex  entropy 
function.  On  the  other  hand,  Conley  and  Saoller  [1]  have  discovered  puzzling  exaaples  of 
constant  positive  definite  aatrlces  D  which  are  inadalaslble,  so  that  for  soam  family  of 
shock  waves  the  systea  of  ODK's  (1.6)  fails  to  adalt  a  trajectory  satisfying  (1.7). 
Moreover,  these  exaaples  exist  in  a  simple  context,  that  arising  when  (1.1)  represents  the 
equations  of  isentroplc  gas  dynamics 


(1.8) 


T  -  V 
t  X 


Vfc  ♦  p(T) 


where  on*  assumes  that  p’ (T)  <  0,  p*(T)  >  0.  The  system  (1.8)  is  also  called  the 

p-ayatea. 

The  natural  requirement  of  parabolicity  imposed  by  Galfand  on  (1.5)  neana  that  tha 
eigenvalue*  of  D( uQ )  alwaya  have  poaltlve  real  parta.  Equivalently,  if  x^(C), 

J  “  1,...,a,  are  the  eigenvalue*  of  the  aatrix  syabol 

(1.9)  PIC)  -  -iCA(ufl)  -  C2vd(u0) 

obtained  froa  linearising  (1.S)  at  a  constant  state  u^,  then  for  aoae  4  >  0  depending 
on  Ug  these  eigenvalues  satisfy 

(1.10)  Re  x^C)  «  -6|C|2  for  ICI  large.  )  -  1 . a  . 

Our  central  objective  in  this  paper  is  to  introduce  a  very  natural  algebraic 
requirement  on  the  linearised  ayatea  froa  (1.5)  beyond  the  condition  (1.10)  -  a  condition 
on  the  viscosity  aatrix  D(u0)  we  cell  strict  stability.  Given  the  concept  of  strict 
stability,  the  bulk  of  this  paper  (sections  2-4)  Is  devoted  to  two  goelsi  1)  To  explain 
the  exaaples  of  inadmissibility  in  [1],  and  link  the  aachanlaas  of  inadmissibility  with 
quantitative  violation  of  the  strict  stability  condition!  2)  TO  elucidate  the  close 
relationship  between  strict  stability  and  necessary  and  sufficient  conditions  for 
admissibility. 

The  requirement  of  strict  stability  ia  motivated  by  tha  following  considerations.  The 
main  Interest  in  the  problea  of  finding  viscous  shock  profiles  as  described  in  (1.5)— (1.7) 
is  to  investigate  in  a  special  case  the  limit  as  v  ♦  0  of  solutions  of  (1.5).  The 
constant  states,  Uq,  are  special  solutions  of  these  diffusion  equations.  Linearization 
of  (1.5)  around  the  constant  state  uQ  yields  the  linear  equation. 


(1.11) 


vt  +  A(u0,vx  *  vD(VTxx'  v  51  0 


v  (x,0)  -  v0(x) 


so  that  a  natural  minimum  requirement  for  any  viscosity  aatrix  in  (1.5)  is  that  the 


-3- 


This  Is  true  (or  all 


solution  vW  of  (1.11)  converges  to  v°  for  sny  initial  data  Vg. 

Vg(x)  €  L2(R)  if  and  only  if  for  any  T  >  0,  there  is  a  fixed  constant  C(T)  so  that 
with  Sv(t)vfl  -  vV, 

(1.12)  Bax  ISV(t)v  I  <  C(I)I»  I  . 

0<w<1  L  L 

0<t<T 

That  ia,  the  initial  value  problem  (1.11)  is  required  to  be  unlforely  well  posed  in  L2  as 

v  ♦  0.  At  any  given  value  Ug,  there  is  a  sat  of  a  x  b  Matrices  D(ug),  the  uniformly 

stable  Matrices.  S(Ug),  guaranteeing  (1.12).  However,  this  set  of  Matrices,  S(ug),  is 

a  bit  too  large  since  the  boundary  of  S(uQ),  3s (uQ),  Includes  0=0  as  well  as  n  x  n 

Matrices  for  which  the  solutions  of  (1.11)  have  a  purely  dispersive  character  (see  section 

4).  The  set  of  strictly  stable  viscosity  Matrices  at  the  point  uQ  is  the  interior  of  the 

eet  S(ug).  The  strictly  stable  viscosity  Matrices  at  tig  adait  the  following  algebraic 

characterisation  (see  section  2  and  [13]  for  further  results)! 

A  viscosity  Matrix  ia  strictly  stable  if 
and  only  if  there  exists  a  S  >  0  so  that 

(1.13)  the  eigenvalues  k^(5),  1  <  j  <  b  for 

the  synbol  P(5)  “  -A(Ug)iC  -  ?2D(Ug)  satisfy 
Re  *  (O  <  -5 |C f 2  for  all  [Cl  . 

Looking  back  at  (1.10),  we  see  that  the  condition  of  strict  stability  strengthens  the 
requirement  of  parabollcity  to  an  algebraic  stability  condition  valid  for  all  {  e  »'  and 
not  just  for  sufficiently  high  wave  numbers.  Thus,  one  objective  here  is  to  study  the 
existence  of  viscous  shock  profiles  for  diffusion  matrices  satisfying  (1.13)  for  every 
value  of  Ug  —  we  call  these  the  strictly  stable  viscosity  natrices  in  the  remainder  of 
this  paper.  Our  second  objective  is  to  explain  the  inadmissibility  example*  from  [1]  and 
to  identify  the  concrete  mechanisns  of  lnadaissibility  through  violation  of  the  strict 
stability  conditions  in  (1.13). 

In  section  2,  first  we  study  sobs  of  the  algebraic  implications  of  strict  stability 
regarding  the  linearised  structure  of  the  OCR's  in  (1.6)  at  the  critical  points,  u^, 
for  a  general  a  x  m  system.  For  a  k-shock  solution  of  (1.1)  satisfying  Lax's  entropy 
Inequalities, 


-4- 


V'iU 


(1.14)  >  »  >  *fc(“R)»  Xlc+1<UR*  >  •  >  Vl'V 
we  verify  there  that  for  the  ODE  in  (1.6), 

(1.15)  the  dimension  of  the  unstable  manifold  at  Is  ■  -  k  ♦  1 

the  dimension  of  the  stable  manifold  at  uR  la  k 

The  following  Is  an  Immediate  consequence  of  this  simple  faett 

Corollary  1.  All  the  explicit  inadmissible  viscosity  matrices  for  the  system  In  (1.10) 
constructed  in  [1]  (via  Theorems  3.2  or  S.2  in  (1))  are  not  uniformly  stable  at  either 

"L  or  V 

Thus,  the  natural  requirement  of  strict  stability  discussed  in  (1. 1 1)-( 1. 13)  is  violated  In 
these  explicit  examples  of  insdmlssibllity .  In  section  2  we  also  develop  apparently  weaker 
algebraic  conditions  which  are  equivalent  to  the  strict  stability  condition  In  (1.13)  for 
2x2  systems  —  these  results  are  central  in  our  discussion  of  the  isentroplc  gas 
dynamics  aquations  in  section  4. 

In  section  3  we  study  weak  shock  profiles  for  general  m  x  m  systems.  We  prove  a 
theorem  which  essentially  characterises  matrices  admissible  for  weak  k -shocks.  An 
immediate  consequence  is 

Corrollarv  2.  Assume  that  D(u0)  is  strictly  stable  at  Uq  for  the  m  x  m  system 

(1.1) .  Then  there  is  a  fixed  neighborhood  of  Ug  such  that  for  any  weak  solution  (1.2)  of 

(1.1)  (satisfying  (1.3))  with  and  u^  in  that  neighborhood,  u^  and  can  be 

connected  in  (1.6)  by  a  viscous  shock  profile  satisfying  (1.7)  if  and  only  if  the  weak 
solution  (1.2)  satisfies  Liu's  strict  entropy  condition  (see  section  3). 

Our  proof  of  this  theorem,  based  on  the  center  manifold  theorem,  is  quite  simple  and  was 
motivated  by  the  work  of  Kopell  and  Howard  [8],  which  might  be  applied  in  the  genuinely 
nonlinear  case.  However,  some  new  observations  are  needed  to  handle  the  general  (non- 
genulnely  nonlinear)  case.  So  strictly  stable  viscosity  mstrices  D(u)  are  always 
admissible  for  sufficiently  weak  shocks.  On  the  other  hand,  for  2x2  systems  we  have 
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Corollary  3.  Suppose  ■  -  2.  It  the  2x2  viscosity  matrix  D(u)  is  not  stable  at 


Uq,  necessarily  it  is  Inadmissible  for  all  weak  k-ahocks  with  elthsr  k  -  1  or  k  «  2. 

In  section  4  we  present  new  examples  of  inadmissible  viscosity  matrices  in  the  large 
where  no  connection  of  and  uR  is  possible,  satisfying  (1.7)  or  the  reverse.  In 

these  examples,  violation  of  strict  stability  is  linked  directly  to  Hopf  bifurcation  in  a 
way  which  elucidates  the  mechanism  behind  the  nonconstructive  Theorem  5.3  in  [1] . 

All  the  above  results  might  lead  us  to  guess t 

Strictly  stable  viscosity  matrices,  for  reasonable  2x2 
genuinely  nonlinear  systems  are  always  admissible  for  any 
shock  solutions  of  (1.1)  satisfying  tax's  shock  inequalities. 

In  section  4  we  identify  broad  classes  of  strictly  stable  D(u)  which  are  admissible  for 

all  shocks  of  the  p-system  (1.8).  We  mention  here  same  special  results: 

Corollary  4.  Consider  the  p-system  (1.8),  and  assume  p'(t)  ♦  0  as  t  ♦  p'(t)  ♦  -• 

as  r  +  0. 

1)  Any  constant  strictly  stable  viscosity  matrix  D  for  (1.8)  is  globally 
admissible. 

2)  Given  any  viscosity  matrix  Dfl  strictly  stable  at  a  fixed  point  uQ  e  R2 ,  there 
is  a  smooth  strictly  stable  D(u)  such  that  D(uQ)  -  D0  and  D(u)  is  globally 
admissible. 

Our  proofs  make  use  of  Lyapunov  functions  obtained  from  convex  entropies  in  a  way  suggested 
by  Mock's  work  (12). 

Despite  these  positive  results,  the  conjecture  above  fails.  We  can  construct  a 
(rapidly  varying)  strictly  stable  D(u)  for  the  system  (1.8)  which  is  inadmissible  for  any 
given  shock.  However,  the  mechanism  of  inadmissibility  is  subtle  and  very  different  from 
those  discussed  earlier. 


Also  In  Motion  4  we  present  an  example  with  D(u)  e  3S(u)  foe  all  u  (stable  but 
not  strlotly  stable)  foe  tha  p-syete*  (1.8)  for  which  (1.5)  la  a  dispersive  system  and 
(1.8)  la  a  conservative  system  admitting  a  first  lntagral,  and  no  shock  profllaa. 

Finally/  wa  remark  that  tha  eholca  of  tha  l2-notm  for  dataralnlng  a  class  of  llnaarly 
atabla  viscosity  matrices  Is  not  tha  only  natural  choice.  Other  natural  norma  for  shock 
wave  theory  such  as  tha  L1  or  BV  norms  might  single  out  a  smaller  class  of  stable 
viscosity  matrices  which  perhaps  admit  shock  profllaa  with  more  special  structure. 


{2.  THE  ALGEBRAIC  8TRUCTURE  IMPLIED  BY  STABLE  VISCOSITY  MATRICES 

If  we  apply  Fourier  transforms  and  Plancherel's  theorem  to  (1.13),  we  conclude  that  a 

viscosity  matrix  is  uniformly  stable  if  and  only  if  the  uniform  bound 

(2.1)  max  |exp  ^  P(OI  <  C(T> 

0<t<T 
0<  v 

is  satisfied  5  e  R  where  ((()  ”  -t2D(uQ)  -  iCA(uQ).  In  appendix  A,  we  prove  the 
following  algebraic  characterization  of  the  strictly  stable  viscosity  matrices,  using  the 
Kreiss  matrix  theorem: 

Theorem  2.1. 

The  following  are  equivalent  for  a  viscosity  matrix  D(u)  for  the  m  x  m  system  in  (1.5) 

1 )  D( u)  is  strictly  stable  at  Uq 

2)  The  eigenvalues,  k^(£)<  j  ■  1,...,m  of  P(£)  satisfy 

Re  x^(5)  <  -fiQ|5|2  for  some  fixed  4Q  >  0  and  all  (  (  t  , 

3)  The  following  three  conditions  are  satisfied: 

(i)  The  system  in  (1.5)  is  parabolic,  i.e.  the  eigenvalues  of 
D(uq)  have  positive  real  part 
(il)  t^DTjJUg)  >0,  k  -  1,...,a 

(ill)  The  symbol  P(£)  has  no  purely  imaginary  eigenvalues  for  5  *  0 
Remark.  We  note  here  that  a  simple  sufficient  condition  guaranteeing  strict  stability  is 
the  following:  There  is  a  positive  definite  symmetric  matrix,  E(Uq),  so  that  EA(uq)  is 
symmetric  and  EO(uq)  is  positive  definite  (perhaps  not  symmetric). 

For  the  proof  of  this  fact,  we  compute  that 

Re(x)e*Ee  ♦  «2Re  e*SDe  -  0 

where  (P(5)  -  *)e  “  0,  so  that  criterion  (2)  of  Theorem  2.1  is  satisfied. 

Next  we  use  theorem  2.1  to  determine  the  linearized  structure  of  (1.6)  at  the  critical 
points  Uj,.  ur.  If  ug  is  any  critical  point  of  (1.6),  the  stable  (unstable)  manifold 
W_(M+)  of  (1.6)  at  Uq  is  tangent  to  the  invariant  subspace  of  the  linearization  of  (1.6) 

St  Uq, 


■8- 


(2.2) 


D_1(A(u0)  -  *1)  »  Q( Uq , s ) 
corresponding  to  eigenvalues  with  negative  (positive)  real  parts,  and 

(2.3)  dim  M_(u0)  3  number  of  eigenvalues  of  Q(uq,s) 

with  negative  real  parts 
dim  M+!u0)  -  number  of  eigenvalues  of  Q(Uq.s) 
with  positive  real  parts 


We  have  the  following  general  fact: 

Theorem  2.2. 

Suppose  D(Uq)  is  a  strictly  stable  m  »  m  viscosity  matrix  for  (1.1).  For  any  k-ahock 
satisfying  Lax's  entropy  inequalities. 

W  >  8  >  VV 


WV  >  8  >  Vi,Ul' 


it  follows  that 


(2.4)  dim  M_(uR)  -  k 

dim  M+(Uji)  -  m  -  k  +  1 

Before  proving  Theorem  2.2,  we  remark  that  Corollary  1  of  the  Introduction  follows 
immediately  in  the  following  fashion:  For  shocks  moving  with  positive  wave  speed  for  the 
p-system  (1.8),  k  »  2,  so  that  applying  Theorem  2.2  we  have 

dim  M_(ur)  “  2,  dim  M+(Uj<)  *  1 

for  strictly  stable  viscosity  matrices.  On  the  other  hand,  all  of  the  inadmissibility 
criteria  in  Theorem  3.2  and  Theorem  5.2  of  (1)  imply  that  necessarily 

dim  M_(uR)  »  1,  dim  M+(uL)  "  2  . 

Thus  the  Inadmissible  viscosity  matrices  constructed  through  these  criteria  cannot  be 
strictly  stable.  Similar  remarks  apply  for  shocks  in  the  p-system  with  s  <  0,  and  for 
applications  of  the  results  in  [1]  for  general  2x2  systems. 

Proof  of  2.2. 

We  consider  the  m  +  1  open  intervals, 

x0  * 

rj  “  aj(u0)'Xj-M<U0,)'  3  ”  ’ . .  -  '•  \  “  ‘  W'-*  ’ 
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M«  nt 


»*  -  number  of  eigenvalues  of  Q(ug,a)  with 
positive  real  parts 

s'  •  number  of  eigenvalues  of  Q(Ug,s)  with 
negative  real  parts 

and  begin  the  proof  of  Theorem  2.2  with  the  claim, 

(2.S)  s\  a*  are  constant  as  s  varies  over  each  1^ 

and  »*  *  s“  •  m 


provided  that  D(ug)  is  strictly  stable.  First,  from  (1.14),  zero  is  never  an  eigenvalue 
on  any  I  y  furthermore, 

(2.6)  Q(u g,s)  never  has  any  non-zero  purely 

imaginary  eigenvalues,  it  *  0,  for  any  s 


because 


det(Q  -  it)  ■  0  <»■>  det(P(T)  -  k)  »  0 

with  *  m  -iar  so  that  k  is  on  the  imaginary  axis)  this  contradicts  iii)  of  3)  in 
Theorem  2.1.  Since  the  eigenvalues  of  Q  are  continuous  functions  of  s  and  cannot  cross 
the  imaginary  axis  on  1^,  we  deduce  the  claim  in  (2.5).  The  proof  of  Theorem  2.2  is 
finished  once  we  establish  that 

(2.7)  s^VU  “  “  -  J  "  °.t,...,m 


for  any  Ug  (apply  the  Information  from  Lax's  shock  Inequalities  and  (2.6)  with 

Ug  «  u^).  As  s  ♦  «•,  Q(Ug,s)  ♦  -sD  1  and  3)  i)  of  Theorem  2.1  together  with  (2.5) 

guarantees 


(2.8) 


0 


From  (2.6),  we  see  that  an  eigenvalue  of  Q(ug,s)  can  cross  the  imaginary  axis  only 


through  sero  and  at  the  special  points,  s  -  l^tug),  k  *  1,...,m. 

For  |s  -  A^l  <  e,  q(uq,s)  has  a  simple 
eigenvalue  Tfc(s),  with  “  0,  and 


-(W 


-1 


k  ■  1 , . . . ,m 


We  claim  that 
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(2.9) 


Fro*  (2.8),  (2.9),  and  3)  ii)  of  'Rieorea  2.1  we  deduce  (2.7)  by  induction  on  j.  It 
remains  for  ua  to  establish  (2.9).  Now,  -A^I  +  A  has  simple  eigenvalues)  therefore, 

-TO  *  (-si  +  A)  has  a  unique  saoothly  varying  slsple  eigenvalue,  8(t,s),  with 
9(0, A^)  ■  0  for  I t |  ♦  |s  -  i  <  e.  Me  let  r(i,s)  denote  the  corresponding  right 
eigenvector  with  r(Q, A^l  -  r^  ( uQ )  satisfying 

(-TO  ♦  (-si  ♦  A) )  r  -  8(1, sir 

By  differentiating  this  formula  with  respect  to  i,s,  evaluating  at  s  -  A^,  t  »  0,  and 
taking  the  inner  product  with  ik(uQ)  we  obtain 


(2.10) 


29 


39, 


q  o  ^ |  pv 

3s  .  ,  '  3i  , 

(0,  A.  )  (0, 


V 


-‘VV  *  0 


39 


where  we  have  applied  3)  ii)  of  Theoree  2.1  to  — .  Froe  the  implicit  function  theorem  and 

(2.10) ,  we  deduce  (2.9)  directly. 

The  Special  Structure  of  stable  2«2  Systems 

First,  we  prove  that  the  sufficient  condition  mentioned  below  Theorem  2.1  is  also 
necessary  when  a  *  2.  In  fact,  we  have 

Proposition  2.1.  The  following  are  necessary  and  sufficient  conditions  that  D(u)  be 
strictly  stable  at  Ug  when  a  »  2. 

1)  D'1  is  strictly  stable  at  Ug. 

2)  There  exists  a  smoothly  varying  positive  definite  symmetric  matrix,  B(u), 
defined  near  uQ  with  SA(u)  symmetric  and  ED(u)  positive  definite. 

3)  (i)  (t^.Dr^) (Ug)  >  0  for  k  -  1  and  2,  and 

(ii)  det  D(u0)  >  0. 

Given  the  matrix  A(u)  with  distinct  eigenvalues,  we  let  R(u)  denote  a  fixed 
saoothly  varying  right  eigenvector  matrix  with  L(u>  the  corresponding  left  eigenvector 
matrix  (LR  «  I,LAR  diagonal).  The  sufficient  condition  for  strict  stability  in  the 
remark  below  Theorem  2.1  is  clearly  satisfied  provided 

(2.11)  LOR  is  positive  definite  . 

But  R(u)  is  not  uniquely  determined  by  A(u).  Proposition  2.1  follows  easily  from  the 
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following  result)  which  says  that  (2.1t)  is  necessary  for  strict  stability  in  2*2 

systeaa  t 

Proposition  2.2.  Assume  a  -  2.  Let  D(u)  be  a  smooth  2x2  matrix.  Then  D(u)  is 
strictly  stable  at  each  point  if  and  only  if  there  exists  a  smooth  positive  diagonal 
matrix  s(u),  such  that  if  R  “  R S,  L  «  R  ' ,  then  LDR(u)  is  positive  definite  for 


We  leave  it  as  an  exercise  for  the  reader  to  Chech  Proposition  2.2  implies  Proposition 
2.1.  The  only  if  part  of  Proposition  2.2  is  the  key  fact  and  this  follows  from  the  lemma 
below. 

Lsmaa.  Suppose  A  -  (*  satisfies  a,d  >0,  ad  -  be  >  0.  Por  -1  <  a  <  1,  let 

diag(/1  -  a,  /I  +  a).  Then  there  exist  a_  and  a+,  -1  <  a_  <  a+  <  1,  depending 

smoothly  on  A,  such  that  for  any  a  with  a  <  a  <  a  ,  the  scaled  matrix  S  1 AS  is 

•  ♦  a  a 

positive  definite. 

Proof.  (“  is  positive  definite  if  and  only  if  a, 6  >  0  and  ( B  *  Y)2  -  4ai  <  0.  Now 

Sa1ASa  -  ic*x  £*)#  *  “  For  -1  <  a  <  1,  this  is  positive  definite  exactly 


Q(o)  =  (b(1  +  o)  +  c(  1  -  a))2  -  4ad(  1  -  o2)  <  0 
To  prove  the  lemma,  we  shall  show  that  Q  <  0  between  distinct  roots  a_  and  a+  of 


Q(a)  with  -1  <  a_  <  a+  <  1.  He  compute 


Q(-1)  -  4c  >  0 


Q(  1 )  «  4b  >  0 


Q(0)  -  (b  +  c)  -  4ad  Q'{0)  .  2(b  ♦  c)(b  -  c) 

J  fi"  *  4ad  +  <b  -  c)2  -  4 (ad  -  be)  +  (b  +  c)2  >  0 


The  minimum  of  Q(a)  is  thus  attained  at  •  -fi’(0)/Q".  Either  |b  -  c|  <  |b  ♦  c|  or 
vice  versa,  so  ( aQ |  <  1.  The  discriminant  of  fl  is 

7  Q’<0)2  -  |  Q"Q(0)  «  <b  +  c) 2(b  -  c)2  +  ( 4ad  -  (b  +  c)2)(4ad  ♦  (b  -  c)2) 


4ad(4ad  -  4 be)  >  0 


v '' v- 


r 


So  Q(a)  has  distinct  rsal  roots  in  the  interval  [-1,1],  which  are  then  smooth  functions 
of  the  coefficients. 

The  technical  result  below,  which  characterizes  the  variation  possible  when  a  vector 
field  is  multiplied  by  an  arbitrary  smooth  2  *  2  strictly  stable  matrix,  will  be  useful 
in  our  investigation  of  inadmissibility  for  strong  shocks  in  section  4.  As  usual,  assume 
that  R(u)  ”  (rj,r2)  is  a  smooth  matrix  of  right  eigenvectors  of  Mu). 

Proposition  2.3.  Suppose  nonzero  vectors  w  and  w  in  R2  satisfy  the 
condi tion 

(Q)  There  exists  a  positive  diagonal  matrix  S  so  that  w  Sw  >  0. 

(Roughly,  condition  (Q)  means  that  w  and  w  are  not  of  opposite  sign.)  Let 
v  m  wtrt  +  w2r2  *  *[“)*/  v  ”  R(u)w.  Then  there  exists  a  strictly  stable  matrix 
D(u,w,w),  depending  smoothly  on  u,w  and  w,  such  that  5v  -  v.  Conversely  if  D  is 
strictly  stable  and  w  any  nonzero  vector  in  R2,  then  w  and  LDRw  satisfy  condition 
(Q). 

The  proof  is  postponed  until  Appendix  B. 

Remark.  Condition  (Q)  covers  two  cases i 

Ql)  If  w  lies  on  a  coordinate  axis  (either  w,  or  w2  are  zero  so  v  -  rA) ,  then 

w  must  lie  in  an  adjacent  quadrant  (so  v^w  >0). 

Q2)  If  w  lies  strictly  in  some  quadrant  Q  of  R2  (w1,w2  nonzero)  then  w  must 

not  lie  in  the  closed  quadrant  opposite  (w  f  -Q). 

Re—rk.  If  w,w  satisfy  condition  (Q),  so  do  w,tw  ♦  (1  -  t)w,  for  any  t, 

0  *  t  4  1,  and  so  do  w,  Sw  for  any  positive  diagonal  S. 
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§3.  ADMISSIBILITY  XH  GENERAL.  FOR  g«  k- SHOCKS 

Ht  begin  this  Motion  by  defining  Liu's  strict  condition  (E).  He  consider  the 
structure  of  the  Hugonlot  set  of  pairs  of  vectors  (u^u^)  satisfying 
(3.1)  tiMj)  -  flu^)  -  »(uL  -  u^)  -  0 

for  soae  wave  speed  s.  Fixing  u^,  the  locsl  structure  of  the  set  of  states  uR 
satisfying  (3.1)  is  well-known  (see  [3]  and  [10]).  In  sosm  neighborhood  of  this  set 

consists  of  m  curves.  u*(p),  k  “  1,..,»  passing  through  u^  with  corresponding  shock 
speeds  sk(p)  for  k«  1,...,e  satisfying 


Theorem  3.1.  Fix  Ug  e  if*  and  k,  1  <  k  <  a.  kssume  X^tu)  la  not  linearly  degenerate 
In  any  neighborhood  of  Ug.  Assume  that  D(ug)  aatlafloa  the  nondegenaracy  conditions: 

(1)  D(Ug)  la  nonsingular 
(11)  lkDrk<u0)  *  0 

(ill)  t-e20  ♦  lUA  -  X^HtUg)  la  nonalngular  for  all  roal  C  *  0. 

Than  tho  following  are  equivalent: 

1)  W%)  >  0  IWV  <  01 

2)  0  1*  locally  {in)afeissible  for  k-ahock*  in  a  neighborhood  of  Ug.  That  la, 

there  axlata  4  >  0  ao  that  for  any  u^  and  in  Bg(Ug)  •  {u  1  |u  -  uQ|  <  4} 

It 

aatlafing  the  Rankine-Rugoniot  relation*  for  some  speed  a  »  a  (pR),  then  a  shock  profile 
lying  in  Bg(Ug)  exists  connecting  u^  to  uR  tuR  to  u^l  If  and  only  if  Liu's  strict 
entropy  condition  s(B)  is  satisfied.  In  any  case,  at  aost  one  trajectory  u( £)  of  (1.6) 
connecting  uR  and  u^  exists  which  remains  in  Bg(Ug)  for  all  £  real. 

Remark i  Before  diecussing  the  proof  of  Theorem  3.1,  we  Indicate  the  fashion  in  which 
Corollary  3  of  the  Introduction  is  an  immediate  consequence  of  Proposition  2.1  and  this 
theorea.  From  3)  of  Proposition  2.1,  if  a  parabolic  viscosity  aatrix  D(u)  is  not 
uniformly  stable  at  Ug,  then  necessarily, 

either  IjOr^Ug)  <  0 
or  tjDr^Ujj)  <  0 

(or  both  occur).  Kith  the  above  inequalities,  we  apply  the  inadmissibility  criterion  from 
Theorem  3.1  to  either  the  one  waves  or  two  waves  to  deduce  Corollary  3. 

Theorem  3.1  is  proved  in  two  steps.  First,  for  all  u^  near  Ug  and  s  near 
X^Ug)  we  reduce  the  connection  problesi  for  the  system  (1.6)  to  that  for  a  scalar  ODE 
locally,  by  employing  the  center  manifold  theorea  with  the  nondegeneracy  conditions 
1-iil).  That  is,  a  curve  is  constructed,  locally  Invariant  for  (1.6),  which  contains  all 
the  critical  points  of  (1.6),  for  any  u^,  in  a  fixed  neighborhood  of  Ug.  In  the  second 
step,  this  one-dimensional  flow  la  analysed:  Critical  points  on  the  Invariant  curve  are 
points  \lo)  on  the  Hugoniot  curve  for  u^  having  s  (p)  *  s.  The  stability  of  the  rest 
point  u^  in  the  flow  is  determined  by  the  sign  of  i^Dr^  ( u^ )  ( X^  ( u^  >  -  *)•  For  a  shock 
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satisfying  s(E),  sk(0)  -  X  (u  )  >  a.  Daganarata  cases  are  treated  by  continuity,  using 

1C  L 

the  center  oumifold. 

Step  1.  Extend  the  systea  (1.6)  by  introducing  the  parameters  v  *  uL  and  s  as 
additional  variables)  then  (1.6)  may  be  written 

Uj  »  D~1(u)tf(u)  -  f(v>  -  s(u  -  v)J 
(3.3)  v^  »  0 

*5  “  0 

Our  analysis  will  be  based  on  the  construction  of  a  center  Manifold  for  (3.3)  at  the 
critical  point  (u,v,s)  “  (u0,u0,Xk(u0) ) .  The  center  manifold  theorem  (Kelley,  [6], 
Theorem  3)  says: 

Theorem.  Suppose  that  a  systea  of  ordinary  differential  equations  aay  be  written  as 

x'  -  Ax  +  X(x,y,  s) 
y'  -  By  ♦  Y(x,y,  z) 
s'  -  Cz  ♦  Z(x,y,z) 

where  A,B,  and  C  are  constant  square  matrices  whose  eigenvalues  have  positive,  zero, 
and  negative  real  parts,  respectively,  and  X,Y,  and  Z  are  Cr(r  >  2)  and  vanish  along 
with  their  first  derivatives  at  (x,y,z)  «  0. 

Then  there  exists  a  locally  invariant  manifold  for  this  systea, 

H*  ■  { ( x,y, z ) |  |y|  <  5,  x  -  u*(y),  z  -  w»(y)} 
where  u*  and  w*  are  Cr  functions  defined  for  |yf  <  5  for  some  6  sufficiently 
small,  and  vanishing  with  their  first  derivatives  at  y  ■  0. 

The  center  manifold  need  not  be  unique,  but  the  following  uniqueness  property  for 
trajectories  does  holds  The  center  manifold  (parametrized  by  y)  may  be  taken  to  be  the 
intersection  of  a  center-stable  manifold  (parametrized  by  y  and  z)  and  a  center- 
unstable  manifold  (parametrized  by  y  and  x).  Then  any  trajectory  which  lies  in  a  small 
neighborhood  Bg(0)  for  all  time  must  lie  on  this  center  manifold.  This  property  follows 
from  this  fact,  stated  in  Kelley  [7] s  If  a  trajectory  starts  in  a  small  neighborhood 
Bg(0)  at  a  point  not  on  the  center-stable  manifold,  then  it  must  leave  B ^ ( 0 )  at  some 
positive  time. 
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Let  ua  now  apply  the  cantar  manifold  theorea.  Without  loss  of  generality,  wa  say 
asauaa  u0  *  0,  X^(Ug)  -  0.  For  convenience,  wa  introduce  w  »  u  -  v  and  the  vector 
W  »  (w,v,a)t  in  K2**1.  Wa  write  (3.3)  In  tha  fora 
(3.4)  W?  -  T(W) 

To  apply  the  canter  manifold  theorea,  it  aufficee  to  deacribe  two  invariant  aubspecea 
for  tha  linearisation  dT  at  tha  critical  point  0t  algebraic  alganapacaa  corresponding 
to  groups  of  aiganvaluaa  with  aero  and  nonsero  real  parts,  respectively.  In  block  fora  on 
W*  x  R*  x  R,  wa  calculate 

r  D-1A(0>  0  0  1 

dT(0 )  -  0  0  0 

L  o  o  o  J 

Tha  characteristic  equation  for  dT(0)  aay  be  written 

X^’detU  -  -  0 

Wa  claia  that  the  nondegeneracy  conditions  iaply  that  the  algebraic  elgenspaoe  for 
eigenvalues  with  zero  real  part  is  siaply  the  kernel  of  dT(0).  Indeed,  condition  iii) 
aaana  that  dT(0)  has  no  nonsero  imaginary  eigenvalues. 

Claia i  The  algebraic  eigenapace  for  dT(Q)  for  the  eigenvalue  0  is  equal  to 
ker  <tt( 0)  if  and  only  if  \Drk(0)  *  0. 

Proof.  This  eigenapace  is  larger  than  ker  dT(0)  if  and  only  if  range  dT(0)  n  ker  dr(0) 
is  nontrivial.  This  occurs  when  there  exists  v  in  ^  such  that  D  *A*v  »  r^,  or 
Drfc  «  A*v,  whence  l^Or^  »  0.  If  t^Dr^  SO,  no  such  v  exists. 

Thus,  let  Y  -  ker  dT(Q),  X  »  range  dT(O).  X  and  Y  are  complementary  invariant 
subapaces  for  dT(0)  comprising  the  algebraic  eigenspaces  for  eigenvalues  with  nonsero  and 
sero  real  parts,  respectively.  Further  explicit  decomposition  of  W2**1  is  unnecessary. 
Applying  the  canter  manifold  theorem,  we  haves 

Proposition  3.2.  Assusm  that  D  satisfies  the  nondegeneracy  conditions  at  ug  •  0  with 
lk<u0>  »  0.  Then  there  exists  4  >  0  and  a  Cr  function  g  i  Y  ♦  X  defined  on 
Bfi(0)  n  Y  so  that 
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1)  m*  -  {x  *  y  e  *-'|x  -  g(y)}  la  a  locally  invariant  manifold  for  aquation  (3.4). 

2)  q(0)  -  0  and  dg<0)  -  0.  Thua  H»  ia  tangant  to  Y  at  0. 

3)  any  trajactory  of  (3.4)  which  liaa  in  1^(0)  for  all  5  liea  in  M*.  In 

particular,  critical  pointa  in  B^{0)  lia  in  M*. 

He  now  deacrlba  how  thia  canter  manifold  reducer  the  connection  problem  for  the  ayatem 
(1.6)  to  one  dl mention.  Obaarve  that  Y  -  ker  dT(0)  la  apannad  by  the  m  +  2  vectora 
(rk,0,0),  (0,0,1),  and  j  -  1,...,m. 

Proposition  3.3.  Xaaume  that  DfUg)  aatiafiaa  the  nondegeneracy  condition a  i-iil).  Than 
there  axiata  6  >  0  eo  that  if  |uL  -  Ogl  +  Is  -  \<u0>l  <  4»  there  la  a  locally 
invariant  curve  u(n,u^,a)  for  the  ayatam  (1.6)  containing  u^  and  any  point  u^  in 
B ^(Ug)  aatiafying  the  Rankine-Hugoniot  relationa  f(uR)  -  ffu^)  -  s(uR  -  u^)  “  0. 

Proof.  Define  a  line  in  Y  parametrized  by  y(n)  -  ( rr^.h^.a).  The  curve 
W(t\)  -  y(n)  *  g(y(n))  liea  in  H»  while  |y(  n)  1  <  S.  Since  g  mapa  into  range 

dT(0),  we  may  write  g(y(n)>  -  (g( n.t^.a) ,0,0 ) .  Under  equation  (3.4),  the  v  and  s 

components  of  W  -  (w,v,a)  remain  invariant.  Hence  the  curve  W(n>  la  the  intersection 
of  two  locally  invariant  manifolds,  ao  ia  locally  invariant.  Returning  to  the  (u,v,a) 
coordinates  of  (3.3),  we  obtain  an  Invariant  curve  for  (1.6)  parametrized  by 

u(n,uL,s)  »  ♦  rrk(u„)  *  g(n.uL,s> 

so  long  as 

ly(n)l  -  |nrkl  +  lu^  -  uQ  |  ♦  Is  -  <  s  • 

If  u„  ia  in  Bs(u0)  and  f(uR)  -  f (u^)  -  s^  -  u^)  -  0,  then  the  point  <u„  -  u^u^s) 
ia  a  critical  point  of  T(w),  so  lies  in  M*,  and  therefore  uR  -  u,^  -  +  g(nR»uL,e) 

for  some  *L. 

The  flow  on  the  invariant  curve  u(n,u, ,s)  ia  now  determined  by  a  scalar  OOB  for 

U 

nU>. 

(3.5)  n5  -  rtn,^,!) 
where  P  ia  Cr,  determined  by  the  equation 

(3.6)  D(u)u^F( n,u^,s)  -  f(u)  -  ff^)  -  a(u  -  u^) 

where  u  •  u(n,u,,s).  Prom  the  uniqueness  property  (3)  of  Proposition  3.2,  two  critical 
L 
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points  in  B4(ufl)  are  connected,  left  to  right,  by  •  trajectory  in  Bj(Ug)  if  and 

only  if  -  0  and  are  connected,  left  to  right,  by  a  trajectory  of  (3.5). 

Step  2.  We  proceed  to  enalyce  the  flow  (3.5).  Two  critical  pointa  ■  0  and  are 

connected,  left  to  right,  by  a  trajectory  of  (3.5)  if  and  only  if  agn  F(n,u^,s)  -  egn 
for  n  between  0  and  n^.  Prca  (3.6)  we  obtain 

P( n.Uj^.a) “  iytt^JIftu)  -  f!^)  -  a(u  -  t^)) 

(3.7) 

V°'V,)  WDtVun  “  l\l\)  -  •>Wun<0'V,> 

we  assume  S  is  so  saall  that  for  each  in  B^tUg)  the  Hugoniot  curves 

-  |r 

u  (p,u^)  in  B^(uq)  are  as  described  at  the  beginning  of  this  section.  The  invariant 
curve  uln.u^s)  intersects  the  Hugoniot  curve  u  (p.u^*  just  when  n  is  a  critical 
point  of  (3.5).  We  define  a  correspondence  between  n  and  p  (given  u^  and  s)  by 
P(n)  -  tJt(uL)(u(n,uL,s)  -  uL),  so  Pn  -  *k<“L>“n 
Tsssis  3.4.  If  6  is  sufficiently  null,  then  if  |u^  ”  “q  |  +  |s  -  <  5,  we  have 

1)  sgn  J^tu^DluJu^  -  sgn  f^Dr^Ug)  and  sgn  -  1  in  B5(uQ).  So  p 

increases  with  n. 

2)  F(n,u.s)  »  0  if  and  only  if  sk(p(n))  “  s  or  n  «  0. 

L 

3)  For  all  n  between  0  and  it. 


(3.8) 


agn  F(n,uL,s)  sgn  I^Or^lUg)  -  sgn  n(s  (p(n))  -  s) 


provided  s  (p(n))  -  a  is  of  one  sign  between  0  and  r^. 

Using  (3.8)  we  nay  complete  the  proof  of  Theoreu  (3.1).  Assume  u^  and 

-1. 

u^  •  u  (p^u^)  satisfy  the  RanKine-Hugoniot  relations  with  s  »  s  (pR),  and  assume  Liu’s 

strict  entropy  condition  s(B)  holds.  Then  uR  «  u ( , u^ , s )  for  some  and 

p_  “  t  (u  )(u_  -  u,  )  «  p( ti  ) .  Ey  (3.8)  and  condition  s(E), 

R  k  L  R  L  R 

sgn  Ftn.u^.s)  sgn  l^Dr^Ug)  -  sgn 

for  all  n  between  0  and  n^.  So  a  trajectory  of  the  flow  (3.5)  connects  0^-0  and 
r^,  left  to  right,  if  and  only  if  I^Dr^lUg)  >  0. 

If  I^Dr^lUg)  >  0  and  u^  is  as  above,  but  the  entropy  condition  is  not  satisfied, 

V 

then  either  s  (p)  *  s  for  some  p  between  0  and  pR,  whence  a  critical  point 
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w-  V  ~ 


Than 


■aparates  0  and  in  (3.8),  or  also  sNp)  <  •  for  all  p  between  0  and  pR. 

(3.8)  implies  that  a  trajectory  of  (3.5)  connacta  on  tha  laft  to  ”  0  on  tha 

right.  In  althar  casa,  no  trajectory  of  (1.6)  lying  in  B^(uQ)  can  connect  on  tha 

laft  to  Ujj  on  the  right. 

Proof  of  Proposition  3.4.  Part  1)  follow*  froa  continuity  and  tha  fact  that 
ur(0,Uq, X^Uq)  )  -  rk(u0>,  since  dg(0,u0 ,  X^tUg ) )  -  0.  Por  part  2),  if  6  is 
sufficiently  snail  and  n  *  0,  than  F(n,uL»s)  "0  if  and  only  if  utn.u^,*)  lias  on  tha 

Swfc  |( 

k-th  Hugoniot  curve  for  u^,  so  u(n,u^,s)  “  u  (p,uL>  for  some  p,  and  s  (p,u^)  “  s. 
(uL)(Sk(p)  - 

sk(p(n)>  >  s  for  n  between  0  and  >  0,  and  ikDrk(uQ)  >  0  (remaining  cases  are 
sisdlar).  First,  X^u^)  “  sk(0)  >  s.  Then  X^tu^)  >  ■  for  any  s  <  a,  so 

>  0  **  close  to  s,  then  p<—  n^u^s)  <  P<r\)#u^vs),  so 

«.  1  U  M.  aw 

F(n,uL,s)  >  0  for  n  between  0  and  j  n^.  (Sines  s  (ptn.u^,*))  >  s,  it  cannot 
vanish  by  part  2).  Letting  a  increase  to  s  wa  gat  F(n,u  ,s)  >  0  for  n  between  0 

4 4 

and  n0-  (Again,  F(n,u^,s)  cannot  vanish  for  n  between  0  and  by  part  2).) 


u^)  “  p(n).  We  shall  establish  part  3)  in  the  case  that 


But  then  p  “ 


Besides  tha  baaic  condition*  p*  (t)  <  0,  p*(t>  >  0,  for  the  p-system  in  (1.8),  In 


thia  section  we  assume  additionally  that 


(4.1) 


-p'(T)  ♦  •  aa  T  0 
-p'<t)  *  0  aa  t  ♦  - 


Thaaa  conditiona  simplify  many  of  tha  statements  below.  Suitable  aodificationa  of  these 
results  when  (4.1)  ia  not  aatiafied  we  leave  for  tha  interested  reader  to  verify,  tilth 
»  (t^.v^),  Ujj  -  ( tr,vr),  the  Bugoniot  relations  from  (1.2)  imply  that 


(4.2) 


-(vL  -  vR)  -  a(tL  -  tr). 


ft 


_p(vj  •  PCtR) 


r  »  *•  '  r  » 


Vi 

) 


p(Tt)  -  p(tR)  -  a(vL  -  vR) 


The  back  shocks  (1-shocks)  are  those  waves  moving  with  speed  a  <  0  and  aatiafylng  the 
entropy  inequality 

(4.3)  "C<V  *  *  *  ~C(  V 

while  the  front  Shocks  (2- shocks)  are  those  waves  saving  with  speed  s  >  0  and  satisfying 
the  entropy  inequality 

(4.4)  c(tL)  >  S  >  c(tR) 

Here  c(T)  »  is  the  Lagrangian  sound  speed,  to  a  consequence  of  Galilean 

invariance,  the  front  and  back  shocks  should  describe  the  same  physics  under  spatial 
reflection.  Indeed,  this  is  the  case  and  in  fact  the  reader  can  easily  verify  that 

(t^.v^),  ( tr,vr)  define  a  front  shock  moving  with 

speed  s  >  0  satisfying  (4.2)  and  (4.4)  if  and  only  if 

(4.5)  (t-»vT),  (t_,v_)  define  a  back  shock  moving  with 

Is  is  K  R 

speed  -s  <  0  satisfying  (4.2)  and  (4.3)  where 
(  t  ,v, )  »  <t_,-v_),  (t_,v_)  -  (  t  ,-v  ) 

Isis  RK  R  R  Isis 

For  th«  p-syfltftm,  onft  right  eigenvector  matrix  is 

<4*s)  *<*>  ■  (dr)  -c!t)) 

with  corresponding  left  eigenvector  matrix 
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lit)  -  J 


Here  we  study  the  (non)exlstence  of  viscous  profile*  for  parabolic  perturbations  of  (1*8) 
where  the  diffusion  aatrix  0  is  given  by 


I 

(  d11 

d12\ 

d11  + 

(4.8) 

D(X,v>  - 

V  d21 

d22  / 

0 

d11d22 

and  tha  aatrix  entries  are  saoothly  varying  functions  of  (x,v).  Proa  (4.6)— (4.8)  we 
coapute 

(4.9)  2  (LOR)  1 1  -  (cd12  ♦  ♦  d11  +  <>22 

2(LDR)22  -  “(cd12  ♦  c_1d21)  +  dn  ♦  A22 

therefore,  froa  3)  of  Propoaition  2.1,  we  conclude  that  D(t,v)  is  strictly  stable  at 
(t,v)  If  and  only  if 

(4.10)  lc(x)d12  ♦  (c(T))_1d21l  <  d^  ♦  d22 
Adalsslble  and  Inadmissible  Viscosity  Matrices  for  Weak  Shocks 

The  following  result  1*  an  isaadiate  corollary  of  Theorem  3.1  and  the  remark  below  that 
theorem: 

Theorem  4.1.  Consider  an  arbitrary  state  <Tg*vg)  and  assume 

<d11  *  d22)  *  <C<Vd12  *  C<V'1d21l  *  °  *  Th*n 

( 1 )  D  is  adalsslble  for  both  front  and  back  weak  shocks  in  a  neighborhood  of 
(t0,v0)  if  and  only  if  D  is  strictly  stable  at  *Tg'vg'" 

(2)  Assume  D  is  not  strictly  stable  at  (tg,Vg). 

A)  If  cd12  ♦  c”1d21  <  -(d-n  ♦  djj)  at  *Tg'vg)  D  18  inadmissible  for  all 
weak  back  shocks  but  0  is  admissible  for  all  weak  front  shocks  in  a 
neighborhood  of  ( x^ , vQ ) . 

B)  If  cd12  +  c-1d21  >  d1t  +  dj2  at  (Tg»Pg)«  D  is  admissible  for  all  weak 
back  shocks  but  D  is  inadmissible  for  all  weak  front  shocks  in  a 
neighborhood  of  ( Xg , ) . 
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We  aika  the  following  two  remark*  which  are  eaay  consequences  of  Theorem  4. 1 


Remark  1.  If  0  la  a  conatant  diffusion  matrix,  than  D  is  admissible  for  all  weak 
shocks  in  a  small  neighborhood  of  all  points  (tq'vq)  and  only  if  D  is  diagonal, 
l.e., 


D  - 


11 


>  0, 


22 


>  0  . 


This  fact  follows  easily  from  Theorem  4.1  since  admissibility  occurs  when 

IcCOd^  +  (c(T))_1d22l  «  d^  +  d22 
and  we  use  (4.1)  with  t  +  0,  t  (  •  to  justify  the  above  remark. 

Remark  2.  There  are  never  any  examples  of  diffusion  matrices  D  for  the  p-system  which 
are  inadmissible  for  all  (front  and  back)  weak  shocks.  The  announced  example  in  the 
introduction  of  (1]  and  described  at  the  end  of  that  paper  is  inadmissible  for  all  front 
shocks  but  admissible  for  weak  back  shocks  at  least  (this  remark  corrects  a  small  error  in 


HI). 

Hopf  Bifurcation  and  Inadmissibility  in  the  Large 

The  examples  of  inadmissibility  described  through  Theorem  3.1  are  not  the  only  ones 
which  occur  in  the  large.  One  of  the  critical  points  can  also  be  encircled  by  a  periodic 
(or  homocllnlc)  orbit  leading  to  oscillatory  behavior  and  preventing  connection  of  the 
critical  points.  This  possibility  was  pointed  out  in  Theorem  5.3  of  [1]  through  a 
nonconstructive  argument.  Here  we  link  the  appearance  of  such  periodic  orbits 
quantitatively  with  violation  of  the  strict  stability  condition.  K  more  subtle  and  quite 
different  example  of  inadmissibility  through  oscillatory  behavior  for  strictly  stable 
diffusion  matrices  is  discussed  later  in  this  section. 

To  bo  specific,  we  assume  that  there  is  a  fixed  point  where  d(Tq/vq>  is 

not  uniformly  stable  (see  (4.10))  and  in  fact, 

(4.11) 

He  consider  front  shocks  with  (tr,vr)  “  *To,wo'  and  the  8hoclt  *P««d,  s,  as  a 
bifurcation  parameter  with  s  >  c(tq).  Prom  (4.2)  and  (4.4)  it  is  easy  to  see  that  given 


d11  *  d22  <  d12C<V  +  d21<C(T0>) 


-1 
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(t  ,v  )  and  ■  >  c(t  )  thara  ia  a  unique  (t  ,v  )  defining  a  front  shock.  Proa  (4.1) 

is  R  A  X*  la 

a  varies  over  (c(tr),»)  ao  that  given  (4.11),  there  ia  a  critical  ahock  wave  speed 
s0  >  c ( )  with 


(4.12) 


^J)  <d11  *  d22> 


C(t0)d12  *  C<V"’d21 


He  compute  that 


(4.13)  tr(D-1('*2 

-c 


-1 


)) 


c 

d*t  D 


[-  f  (dt1  ♦  d22)  4  Cd 


12 


Therefore,  froa  (4.12),  (4.13)  at  the  critical  point  (TR»VR>  “  'To,vo*  and  tor 
la  ”  aQ|  <  5,  D-,(A  -  si)  has  nonzero  coaplex  conjugate  eigenvalues,  1(a),  X(a)  with 


(4.14) 


ReA(a)  >  0  , 
Ra  1(a)  <  0  , 
Re  l(s0)  -  0 


-  S  <  a  <  aQ 
<  ■  <  *0  +  6 


and  also 
(4.15) 


Re  V(„|r|jM 


ao  by  Bopf's  bifurcation  theorem. 


<W 


ia  encircled  by  a  aaall  amplitude  periodic 


orbit  for  either 


•o  -  «• 


<  a  < 


•o  < 


<  a  ♦ 


and  the  corresponding  front  ahock  is  necessarily  inadaissible  for  this  viscosity  aatrix. 


The  saae  phenomenon  at  (t0,vq)  occurs  for  back  ahocka  provided  that  the  strict  stability 
condition  is  violated  through 


d12C(V  *  d21C(V 


<  -(d 


11 


d22> 


rather  than  (4.11). 

By  looking  back  at  Theorems  2.1  and  2.2,  the  reader  can  see  that  the  above  argument 
for  the  p-systaa  illustrates  in  a  special  case,  a  very  general  link  between  violation  of 
strict  stability  and  occurence  of  saall  amplitude  periodic  orbits  bifurcating  froa  critical 
points  for  (1.6)  for  general  a  a  a  systaas  and  guaranteeing  inadmissibility  when  a  «*  2 » 
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►■s.V 


r 


however,  to  keep  our  discussion  brief,  we  do  not  develop  this  here  in  detail  beyond  the 
above  example • 

Admissibility  in  the  Large 

Here  we  exhibit  a  reasonably  wide  class  of  strictly  stable  viscosity  matrices  which 
are  admissible  for  all  shocks  of  the  p-systea.  In  particular,  we  prove  Corollary  4  of  the 
introduction. 

Theorem  4.2,  Suppose  conditions  (4.1)  hold,  and  suppose  D(T,v)  is  a  smooth  strictly 
stabls  viscosity  mstrix  for  the  p-systam  such  that 

a)  d(t,v)  is  constant  exterior  to  a  compact  region  Q  in  the  half  plane  x  >  0. 


/  c(x)2 

b)  For  some  fixed  1,  I 

V  * 


X 

1 


•D(x,v)  is  positive  definite  in  Q. 


Then  d(t,v)  is  admissible  for  all  shocks  of  the  p-system. 

Remark.  For  any  fixed  (x,v),  there  does  exist  X,  |X|  <  c(x),  such  that 
c2  X 

(*  1  )  D( x,v)  is  positive  definite.  To  see  this,  consider  the  matrix  of  left 

eigenvectors  for  the  p-system  defined  by 


/  /TTi 

'l  • 


i 


It  follows  from  Proposition  2.2  that  for  certain  o  in  the  Interval  (-1,1),  the  matrix 
-1  T 

L  DI>  ,  hence  LID,  la  positive  definite.  But 
a  a  a  a 


T 

oi  a 


2  (  ,  )  . 

'co  1  ' 

Corollary  4  of  the  Introduction  follows  immediately  from  the  theorem  and  remark 


above.  For  if  D„  is  strictly  stable  at  a  fixed  ( xQ , vQ ) ,  then  for  some 


( 


c2(x) 


*)  dq  is  positive  definite  at  x  ■  xQ,  so  also  for  x  in  a  small 


neighborhood  0  of  (Xg,Vg).  Let  <|i(t,v)  be  a  function  such  that  t  ■  1  at 
(x0.v0),  ♦  =  0  outside  Q,  and  let  D(x,v)  »  <C0  ♦  (1  -  f)X.  Theorem  4.2  applies. 


yielding  Corollary  4. 
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To  put  our  proof  of  4.2  in  context,  recall  that  an  entropy  for  the  system  of 


conservation  laws  (1.1)  is  a  function  E(u)  such  that  for  some  Q(u)  (the  entropy  flux)  we 

have 


(4.16) 


7E-A(u)  -  Vfi(u) 


Smooth  solutions  of  (1.1)  also  satisfy  E(u)t  +  Q(u)x  -  0.  The  Hessian  of  E(u) 

2 

symmetrizes  (1.1),  l.e.,  7  E*A(u)  is  symmetric  (differentiate  above).  (Compare 

Propoeltln  2.1(2).)  A  2x2  system  such  as  the  p-system  admits  many  entropies.  In  fact, 

2 

for  the  p-system  (4.16)  reduces  to  the  one  equation  Ett  *  c  (x)E  •  0.  He  consider  two 

2  - 

special  solutionsi  E  (t,v)  »  v  /2  ♦  P(r),  where  P'(t)  ■  p(T),  and  E( t,v)  -  tv.  He 
0  2  c2  X 

define  E^  »  EQ  +  XE,  and  note  7  E^  •  ^).  Hypothesis  b)  of  4.2  simply  says  that 

2 

V  E^D  is  positive  definite  in  fi,  and  lsqilies  that  E^  itself  is  convex  in  0. 

Proof  of  4.2.  Plx  a  front  shock  with  u_  -  (t_,v_),  u_  *  (T_,v_)  and  s  satisfying 
“  L  Id  L  a  R  R 

(4.4)  so  that  t  >  t  ,  v  >  v  .  For  the  system 

R  L  L  R 


(4.17) 


D(u)u 


-S(T  -  Tr)  -  (V  -  VR) 


p(T)  -  p ( TR >  -  S(V  -  VR) 


i  V(U) 


we  shall  show  a  trajectory  exists  connecting  u^  on  the  left  to  uR  on  the  right.  From 
Theorem  2.2,  uR  is  a  stable  node  for  this  system  and  la  a  saddle  point. 

He  now  invoke  some  results  of  Conley-Smoller  (1]  and  claim:  If  no  periodic  or 
homocllnlc  orbit  exists  encircling  u^,  then  one  branch  of  the  unstable  manifold  of  u^ 
approaches  uR  as  x  ♦  “ .  The  results  of  [1]  which  are  pertinent  are  a  classification 
theorem  for  flows  in  the  plane  with  two  critical  points,  one  a  saddle,  one  a  node  (Lemma 
4.1),  and  the  existence  of  an  "isolating  disk"  for  the  system  (4.17)  with  a  constant 
diffusion  matrix  (Lemma  5.1),  which  traps  some  branch  of  an  invariant  manifold  of  u^ 
Inside  it  (Lemma  4.2).  In  the  present  situation,  the  isolating  disk  of  Lemma  5.1  may  be 
constructed  to  contain  0  in  its  interior,  since  D(u)  is  constant  exterior  to  ft. 

He  now  Introduce  two  functions 
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(4.18) 


Ao  »  p(T)(s(T  -  tr5  +  w  -  VR)  ~  »(v2/2  ♦  P(t)  -  V(p(TR)  -  SVR) 

A  -  P(T)  -  Tp(TR)  -  S(T  ~  TR)(v  -  VR)  -  (V2/2  -  WR) 

and  define  A^  «  AQ  +  AA.  Thie  A^  coincides  with  the  functional  A  which  appears  in 
Mock  [12]  if  the  entropy  is  taken  to  be  E^.  It  has  the  property  that 

VK\  "  (°  I1’  1^Vlu>  " 

Along  any  trajectory  of  (4.17),  then, 

Ax(u)x  -  V(u)*V2BxD_1V(n) 

2 

So  in  any  region  where  V  S^D  is  positive  definite,  A^  is  Increasing  along  trajectories. 

Mote  that  f roe  the  first  reswrk  of  section  4,  D(t,v)  is  diagonal  exterior  to  9.  It 
follows  that  Aq  Increases  along  trajectories  exterior  to  9,  and  that  A^  Increases 
along  trajectories  in  the  strip  0  <  t  <  where  c(t^)  -  |X|.  (This  strip  contains 

9.)  The  phase  portrait  of  (4.17)  for  D  diagonal  is  shown  in  Pig.  1. 


Figure  1.  Phase  portrait  of  (4.17)  for  diagonal  D 
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He  now  consider  two  cases.  First,  suppose  x  »  x^  Then  it  is  clear  that  no 
periodic  or  hostoclinic  orbit  encircling  uR  can  exist,  for  it  would  have  to  Intersect  the 
vertical  line  t  -  r  In  two  places,  but  rig.  1  shows  that  is  Impossible  for  a  trajectory 
in  the  region  x  >  xR. 

In  the  second  case,  xR  <  x^,  no  periodic  or  homoclinic  orbit  encircling  uR  could 
lie  entirely  to  the  left  of  the  line  x  ■  t  ,  for  A^  increases  on  trajectories  there. 

The  remaining  possibility  is  that  part  of  the  encircling  orbit  lag  to  the  right  of  the  line 
t  •  y  It  can  only  do  so  if  that  part  is  one  oonnected  piece  as  indicated  in  fig.  1. 

Let  u  -  (x^.v  )  be  the  point  of  entry,  u+  “  ( x^ ,v+)  the  point  of  exit.  Then 

AQ(u+)  -  Afl ( u_ )  >  0  >  Axtu+)  -  Ax(u_) 
if  the  encircling  orbit  exists.  Since  1  ■  c(Xj),  this  lilies 

c(  Xj )  |  A(u+)  -  A(u_ )  |  >  A0(u+)  -  AQ(u_) 

He  will  show  this  cannot  hold,  so  the  encircling  orbit  cannot  exist.  Now 

A(u+)  -  A(u_ )  -  (v+  -  v_)(v1  “  (v+  ♦  v_)/2) 
where  V1  “  VR  “  “stx^  -  xR).  with  v2  ”  VR  "  (pt  t1 )  ”  pt  XR)  )/s,  observe  that 
v_  <  v^  <  v+  <  Vj.  Also 

AQ(u+)  -  AgCu)  -  s(v+  -  v_)(v2  -  (v+  +  v_)/2>  . 

Then  v2  -  (v+  ♦  v_)/2  >  |v1  -  (v+  +  v_)/2|  and  s  >  c(XR)  >  c(Xj),  so  the  inequalities 
above  cannot  hold,  concluding  the  proof  of  Theorem  4.2  for  front  shocks.  For  back  shocks, 
replace  v  by  -v,  x  by  -x  in  (4.17),  reducing  the  connection  problem  to  that  for  a 
front  shock. 

An  Example i  Strictly  Stable,  but  Inadmissible  in  the  Large 

He  give  here  a  construction  which  shows  that,  despite  the  positive  results  above,  the 
local  condition  of  strict  stability  is  not  quite  sufficient  for  global  admissibility. 
Proposition  4.3.  Fix  a  front  shock  (xl,v.  ),  (t_,v_)  with  c(  x.  )  >  s  >  c(  x_ )  >  0  There 

'  Lt  Li  K  K  L  R 

eixsta  a  smooth  choice  of  0(x,v),  strictly  stable  at  each  point,  such  that  the  system  of 
ODE's  (4.17)  for  the  p-aystea  admits  a  closed  periodic  orbit  encircling  ( tr, vR) « 
excluding  (x  ,v_  ),  so  no  trajectory  can  connect  the  two.  We  may  choose  DEI  outside 

Is  Is 


an  annular  region  containing  the  periodic  orbit 


Proof,  We  will  sake  iih  of  ths  strictly  stable  matrix  function  D(u,w,w)  of  Proposition 
2.2.  Observe  that  d(u,w,w)  *  I  (see  appendix  b).  Of  course,  u  »  (t,v).  Our  procedure 
is  as  followsi 

a)  He  will  exhibit  a  vector  field  V^tu)  in  on  annular  region  encircling  (tr,vr) 
excluding  ( TL»VL'  which  admits  periodic  orbits,  and  is  also  such  that  the  pair  LV(u), 
tVl(u)  satisfy  condition  (Q)  of  Proposition  2.2  at  each  point. 

b)  Than  v  and  can  be  patched  together  by  a  partition  of  unity i  Take  a 

function  t(u)  which  is  1  on  a  periodic  orbit  and  0  outside  the  annular  region,  and 
let  V  -  ♦  (1  -  *)v. 

c)  By  the  remark  concluding  section  2,  LV  and  LV  satisfy  condition  (Q) 
everywhere.  Simply  take  D-1(u)  »  D(u,LV(u) ,LV(u) ) . 

It  remains  to  perform  step  a) .  Because  u^  “  ( tr,  vR)  is  a  stable  node  for  the  vector 
field  v(u),  the  map  u  ♦  w  given  by 

w  -  -I«V(u) 

is  locally  Invertible  at  w  -  0  by  the  Inverse  function  theorem.  In  the  w-plane  we 
consider  the  family  of  curves  w(8»rfl)  in  polar  coordinates, 

r(«,r0)  -  (I cos  8  -  8Qlp  ♦  lain  8  -  8Q|p)“1/p  rQ 
for  r0  small,  p  end  80  fixed.  1  <  p  <  2,  0  <  8Q  <  w/4  (see  fig.  2) 
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This  curvs  is  ths  boundary  of  a  ball  in  the  p-metric  Iwl^  “  (Iw^l^  +  (w^l**)1'^  rotated 
through  an  angle  8Q.  It  has  the  property  that  for  each  8,  the  pair  of  vectors 
-w(8),  w(9)  satisfy  condition  (Q).  For  rQ  sufficiently  saall,  we  say  then  obtain 
closed  curves  uOjr^)  encircling  Ug  via  the  isomorphism  u  <— >  w  defined  above.  Then 
*<9)  “  i<8>  “  (["  +0C<T>  .  AitJ  +  0<ro))Lu(e,ro) 


-°c(T>^  +  0(r0))^(6,r0) 


There  is  a  saooth  choice  of  positive  diagonal  S(0)  (see  appendix  B)  so  that 
-wT8w(8)  >  CTp  for  all  8,  C  >  0  constant. 

With  s  ■  S'dlagCs  ♦  c,s  -  c),  c  •  c(T(8itQ ))  we  have 

LV(u(8 ) )SLu( 8 )  >  cr*  -  Clr*)  >0  if  rQ  is  saall. 

Therefore  it  suffices  to  take  _  u(8»rQ)  in  step  a). 

The  Admissibility  Criteria  of  conley-Smoller  and  strict  Stablllt 

Here  we  reaark  that  under  very  general  circusstances,  application  of  the  admissibility 
criteria  of  Conley-smoller  frost  [1]  to  one  front  shock  and  one  back  shock  associated  with  a 
given  value  of  T  automatically  forces  ths  diffusion  matrix  to  satisfy  a  stronger 
requirement  than  strict  stability,  ror  simplicity  in  exposition  we  only  state  these 
results  below  for  the  paired  front  and  hack  shocks  described  in  (4.5),  related  through 
reflections!  symmetry,  and  leave  the  obvious  generalizations  to  the  interested  reader.  We 


Proposition  4.4.  Assume  the  admissibility  criterion  of  Theorem  3.2  of  [1J  applies  to  both 
the  front  shock  and  back  shock  described  in  (4.5)  for  a  fixed  shock  speed  s  >  0.  Also 
assume  the  diffusion  matrix  D  satisfies  D(tr,vr)  “  d(tr,-vr)  -  in  particular  any 
constant  D  always  satisfies  this  requirement.  Then  D  is  strictly  stable  at  ( t R,tvR) 
and  satisfies 

-1  C(V 

I  C(  T  )d  +  C(T  )  d  I  *  - —  (<J  +  d  )  <  d  ♦  d 

1  1  «'  12  '  r'  21 1  s  11  22'  11  22 

To  apply  Theorem  3.2  of  [1]  to  a  given  front  shock  connecting  (t ,v  )  to  ( T  ,v  ) 

l  L  R  R 


with  speed  s  >  0  requires  at  (rR,vR), 


d,2s  <  d22 


-V 


•o  that 


c 


21 


<  d 


11 


d12C  *  ^l6'1  <  f  «1,  +  d22> 


Similarly,  connecting  < tR, _VR)  to  ( TL*"VL>  with  wave  apeed  -a  requirea  by  the  au* 

conditiona  that  at  (t_,-v_) 

R  R 

d1SC  +  d21C_1  >  -  f  <d„  *  d22> 

Since  the  entropy  condition  guarantees  c(tr>  <  a,  the  concluaion  of  Propoaition  4.4 
followa. 

We  alao  have 

Propoaition  4.5.  Assume  the  admissibility  criterion  of  Theorem  5.2  of  [11  appiiea  to  both 
the  front  shock  and  back  ahock  deacribed  in  (4.5)  for  a  fixed  ahock  apeed  a  >  0.  hleo 
aeaume  the  diffueion  matrix  D  aatiaflea  DCr^.v^  -  Then  D  ia  strictly 

atable  at  (T^.tv^)  and  aatiaflea 

|c<Vd12  +  C<V‘,d21«  ‘  ^TT  (d11  +  d22> 

I* 

First,  we  apply  the  admissibility  criterion  of  Theorem  5.2  from  [1)  for  the  shock 
moving  with  speed  s  >  0i  this  requirea  that  the  trace  of  the  matrix  in  (4.13)  is  negative 


at  eo  that 


*  C"^d21  *  ld11  +  d22) 
Li 


Similarly,  allying  thia  criterion  to  the  ahock  from  (4.5)  with  speed  -  a  requires  that 

the  trace  of  the  matrix  in  (4.13)  ia  negative  at  (t  ,-v  )  so  that 

L 


Cd12  +  C"’d21  >  <d,1  +  d22» 

la 


and  these  two  inequalities  together  with  the  entropy  condition  c(t^)  >  a  imply  the 
concluaion  of  Proposition  4.5. 

An  Inadmissible  Matrix  on  the  Boundary  of  the  Strictly  Stable  Viscosities 
We  consider  the  explicit  choice  of  the  matrix  D  given  by 


(  ° 

0(T)  -  [ 

U2u> 


1 

0 
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Thia  aatrlx  ia  associated  with  purely  dispersive  wave  propagation  for  (US)  since  LDR(x) 
la  a  skew  sysesetric  Matrix.  Furthermore,  D  ia  on  tha  boundary  of  tha  sot  of  atablo 
viscosity  Matrices  ainca  D( r)  ia  tha  liait  as  e  +  0  of 


and  D£(t)  is  strictly  stable  because  (4.10)  ia  satisfied. 

Fixing  any  shock  (VV'(VV'*'  At  i-  ***y  **  Ch*°k  th*t  “**  function  A0(T,v) 
of  (4.18)  is  constant  along  trajectories  of  (4.17).  Since  A0(TL'VI<>  *  /VTR*TR1'  n° 
connection  is  possible. 
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Appendix  A.  Proof  of  Theorem  2.1 


We  begin  by  developing  some  necessary  criteria  for  a  viscosity  aatrlx  0  to  be 
stable,  l.e.  for  0  e  S(Uq).  From  (2.1)  It  follows  that 
( A.  1 )  max  |e"(C'|  <  C 


0<t<- 

ce* 

tor  som  fixed  constant  C,  where  P{£)  “  -^D  -  i£A.  This  estimate  can  hold  only  if  the 
eigenvalues  of  PH)  have  nonpositive  real  part  for  all  real  (.  Using  this  principle,  we 
can  establish s 

Proposition  (A.1).  Assume  D  is  stable  at  Ug  (D  (  S(Ug)).  Then 
(1)  The  eigenvalues  of  D  have  nonnegative  real  part. 

12)  f^Qr^fUg)  >  0  tor  A  »  1,...,s. 

(3)  For  any  eigenvalue  x^HJ  of  PH), 

M  x  <{)  <  0,  J  -  1 . .  . 

Proof.  From  the  discussion  above,  (3)  is  i mediate.  Define,  for  convenience, 

B(  8)  »  D  Sin  8  +  iA  cos  8  . 

eigenvalues  p^(0)  of  B(0)  are  related  to  eigenvalues  x^H)  of  P (5)  by 
H^( 0) *( tan  8/cos  8)  -  -x^(tan  8)  for  8  a  0,  -w/2  <  8  <  w/2  . 

From  (3),  and  using  continuity, 

(A. 2)  (sgn  0)Re  1^(8)  >  0,  -w/2  <  8  <  w/2,  j  -  . 

Setting  6  -  w/2  we  obtain  (1).  For  (2),  observe  B(0)  -  1A  has  distinct  imaginary 
eigenvalues.  Therefore,  for  small  0  there  exist  msooth  eigenvalues  ( 8)  and 
eigenvectors  1^(8),  with  1^(0)  “  ilj^Ug),  Rfcl0)  "  rh*uo'*  satisfying 

(8<8)  -  u^sniys)  -  0 

Differentiate  and  set  8-0  (B'(0)  -  D) .  Then  dot  with  ( uQ ) .  We  obtain 
(A.3)  tkDrk(u0)  -  6) 

Part  (2)  now  follows  using  (A.2). 

Proof  of  Theorem  2.1.  We  shall  argue  that  (1)  implies  (2)  implies  (3)  implies  (1).  Assume 
that  D  is  in  the  Interior  of  S(Uq).  Then  for  some  >  0,  D  -  4QI  is  stable,  i.e., 
in  Slug).  The  eigenvalues  x^H)  of  PH)  "  -£2(D  -  4QI)  -  HA  then  satisfy 
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— •  tm* 


V  -  V 


?.v  •-  ■&'  ’  Vi.  '• 


Re  SyC)  »  Ra  x^C)  +  «QC2  <  0 

establishing  (2). 

It  is  not  hard  to  show  that  if  sow  part  of  condition  (3)  doaa  not  hold,  then  (2) 
cannot  hold,  by  using  scaling  arguments  as  in  the  proof  of  Proposition  A . 1 .  Since  the 
conditions  in  (3)  are  open  conditions,  to  complete  the  proof  it  remains  only  to  show  that 
if  conditions  31)  -  3111)  are  satisfied,  than  0  is  stable,  i.e,  in  S(Ug),  He  will  make 
use  of  one  part  of  the  Kreiss  Matrix  theorem. 

Theorem  (Kreiss,  1959  [9]).  Let  a  family  of  m  x  m  square  matrices  be  given.  A  necessary 
and  sufficient  condition  that  C3  >  0  exist  so  that 

la1*!  <  C1 

for  all  t  >  0  and  all  A  in  the  family,  is 

(K3)  There  exist  constants  C31  and  C32  and  a  matrix  S(A)  for  each  A  in  the 
family,  with  max(ISI,IS  1 1  )  <  c31<  so  that 


is  upper  triangular,  with 

(A. 4)  Re  IJ  ,<...<  Re  V  <  0 

1  It! 

and 

(A. 5)  Ib^l  <  c32'r*  W j i  for  *  <  j 

Assume  D  satisfies  conditions  31)  -  3111).  He  shall  verify  the  condition  (K3)  for 
the  faatlly  of  matrices  {P(£)K  SR),  or  what  is  the  same,  because  of  a  positive  scaling 
factor,  for  the  bounded  family  {-sin  8  B(8)l-w/2  <  8  <  n/2}. 

We  begin  by  using  condition  ii),  constructing  a  suitable  S  for  8  near  0.  Prom 
the  proof  of  A.1,  the  matrix  B(8)  may  be  diagonalized  by  a  matrix  R(8)  of  right 
eigenvectors  for  |6|  small  (so  S  1  *  R(8),  and  •  0)  and  its  distinct  eigenvalues 

Uk(8)  satisfy  U^(8)  “  t^Dr^tUg)  y  *°  tor  1*1  K  w®  have  -sin  611^(8)  <  0. 
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For  ®0  <  |8|  <  «/2,  we  nay  choose  a  unitary  Mtrix  s(9)  which  puts  -sin  6  b(9) 
into  uppar  triangular  font  satisfying  (A. 4),  by  Schur's  theorem  (or  see  Richtoyer  am 1 
Norton  [14],  p.  77).  The  eigenvalues  u^(9)  of  -sin  9  B(9)  are  continuous  and  never 
touch  the  imaginary  axis  for  8Q  <  |9|  «  »/2  by  conditions  i)  and  iii).  therefore 
| Be  Ujtejl  »  vg  >  0  for  9q  <  |9|  <  r/2.  But  the  off-diagonal  elements  of 
S(-sin  9  B)S  1 (9 )  are  uniformly  bounded,  since  S  is  unitary  and  B(9)  bounded.  So 
(K3)  bolds,  and  theorem  2.1  is  established. 


Appendix  B;  Proof  of  Proposition  2.3. 

Consider  the  cx>nverse  pert  first.  If  0  is  strictly  stable,  there  exists  e  positive 

diagonal  S  such  that  S  ' LDRS  is  positive  definite,  by  2.2.  Then  if  w  is  nonzero, 

1  T — 1  ~  ~-1  T — 2 

0  <  (S  w)  S  LDRS ( S  w)  «  w  S  LDRw  , 

so  w  and  LDRw  satisfy  condition  (Q) . 

Now  assume  w  and  w,  nonzero,  satisfy  (Q).  Then  there  exists  a  rotation  natrix 

0(8)  »  (°?*  q  °in  »),  with  |8|  <  */2,  and  a  constant  c  >  0,  such  that 
sin  w  cos  o 

w  »  c0(8)Sw.  with  LDR  *  C0(8)S,  D  is  strictly  stable,  by  2.1(3),  since  cos  8  >  0.  In 
order  to  show  that  D  may  be  chosen  smoothly,  it  suffices  to  show  that  S,  depending  on 
w,w,  may  be  chosen  smoothly. 

Let  S(t)  -  diag(t,2  -  t).  We  shall  show  that  t  may  be  chosen  as  a  smooth  function 
of  the  angles  8  -  arg  w,  8  «  arg  w  in  the  proper  subdomain  of  the  torus  S1  x  s1.  In 
Figure  3  below,  the  torus  is  divided  into  16  square  patches 
(k  ■J,  (k  ♦  1)  j)  x  (k  j,  (k  ♦  1)  indicating  regions  in  which  w  and  w  lie  in  given 
quadrants.  The  domain  is  the  open,  connected  set  indicated  by  shading. 


Figure  3.  Domain  of  t(8,6)  in  S1  x  S1 


Define  t(8,8)  aa  follows:  If  (9/9)  lies  In  a  square  on  the  Bain  diagonal  (so 

k  »  k,  l.e.,  w  and  w  lie  in  the  same  quadrant)#  define  t(9,8)  «  1,  so  S(t)  “  I. 

Consider  a  particular  patch  0  off  the  main  diagonal,  fi  ■  (o,  *  (j,  »).  Let  w,w  be 

given,  with  <6, 8)  6  Q.  In  order  that  condition  (Q)  above  be  satisfied  with  S  -  S(t), 

we  must  have  t(8,8)  <  t  where  w^Slt  )w  -  0,  i.e.,  w  w  t  +  (2  -  t  )w  w  «  0  or 

C  C  lie  C  2  2 

2T 

- -  »  -tan  8  tan  8  i  T(8,8).  Thus  t  *  - - — 

2  -  t  c  1  +  T 

c 

Simply  take  t(8,8)  »  t(T)  as  any  C°*  function  of  T  on  the  interval  (0,**)  such  that 
t(T)  <  tc  and  t(T)  •  1  for  T  sufficiently  large.  Other  patches  off  the  aain  diagonal 
are  treated  similarly,  so  a  smooth  t(8,8)  may  be  defined  aa  required.  This  concludes  the 
proof  of  Proposition  2.3. 
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We  study  a  natural  class  of  appropriate  viscosity  matrices  for 
strictly  hyperbolic  systems  of  conservation  laws  in  one  space  dimension, 

Uj.  +  f (u)x  «  0,  u  e  R®.  These  matrices  are  admissible  in  the  sense  that 
small  amplitude  shock  wave  solutions  of  the  hyperbolic  system  are  shown  to  be 
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limits  of  smooth  traveling  wave  solutions  of  the  parabolic  system 

+  f(u)x  “  v(Dux)x  as  v  ♦  0  if  D  is  in  this  class.  The  class  is 

determined  by  a  linearized  stability  requirement:  The  Cauchy  problem  for  the 

2 

equation  u^  +  ffug)^  -  MOuxx  should  be  well  posed  in  L  uniformly  in 
v  as  v  ♦  0.  Previous  examples  of  inadmissible  viscosity  matrices  are 
accounted  for  through  violation  of  the  stability  criterion. 
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